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Abstract. The discrete cosine transforms of types V-VIII are generalized to the antisymmetric 
and symmetric multivariate discrete cosine transforms. Four families of discretely and continuously 
orthogonal Chebyshev-like polynomials corresponding to the antisymmetric and symmetric general¬ 
izations of cosine functions are introduced. Each family forms an orthogonal basis of the space of all 
polynomials with respect to some weighted integral. Cubature formulas, which correspond to these 
families of polynomials and which stem from the developed discrete cosine transforms, are derived. 
Examples of three-dimensional interpolation formulas and three-dimensional explicit forms of the 
polynomials are presented. 
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1. Introduction 


This paper aims to complete and extend the study of antisymmetric and symmetric multivariate gener¬ 


alizations of common cosine functions of one variable from 11 . First, the set of four symmetric and four 


antisymmetric discrete Fourier-like transforms from 11 is extended. Then four families of Chebyshev-like 


orthogonal polynomials are introduced and the entire collection of the 16 discrete cosine transforms (DCTs) 
is used to derive the corresponding numerical integration formulas. 


The antisymmetric and symmetric cosine functions of n variables are introduced in 11 as determinants and 


permanents of matrices whose entries are cosine functions of one variable. The lowest-dimensional nontrivial 
case n = 2 is detailed in 9 . These real-valued functions have several remarkable properties such as contin¬ 
uous and discrete orthogonality which lead to continuous and discrete analogues of Fourier transforms. The 
discrete orthogonality relations of these functions are a consequence of ubiquitous DCTs and their Cartesian 
product multidimensional generalizations. There are eight known different types of DCTs based on various 
boundary conditions |3 . Only the first four transforms are generalized to the multidimensional symmetric 
cosine functions [ll . Therefore, the remaining four transforms of types V-VIII need to be developed to obtain 
the full collections of antisymmetric and symmetric cosine transforms with all possible boundary conditions. 
The resulting 16 transforms are then available for similar applications as multidimensional DCTs. Besides 
a straightforward utilization of these transforms to interpolation methods, they may also serve as a starting 
point for Chebyshev-like polynomial analysis. 

The Chebyshev polynomials of one variable are well-known and extensively studied orthogonal polynomials 
connected to efficient methods of numerical integration and approximations. The Chebyshev polynomials are 
of the four basic kinds [5], the first and third kinds related to the cosine functions, the second and fourth 
kinds related to the sine functions. Since both families of antisymmetric and symmetric cosine functions 
are based on the one-dimensional cosine functions, they admit a multidimensional generalization of the one¬ 
dimensional Chebyshev polynomials of the first and third kinds. The results are four families of orthogonal 
polynomials, which can be viewed as the symmetric and antisymmetric Chebyshev-like polynomials of the first 
and third kinds. Several generalizations of Chebyshev polynomials to higher dimensions are known - in fact 
for the two-dimensional case the resulting polynomials become, up to a multiplication by constant, special 


cases of two-variable analogues of Jacobi polynomials 14 16 . In full generality, multidimensional Chebyshev 
polynomials of the first kind are constructed, for example, in 00 . for broader theoretical background see 
also [6,20 . As is the case of the Chebyshev polynomials of one variable, the multivariate Chebyshev-like 
polynomials inherit properties from the generalized cosine functions. The link with these cosine functions 
provides tools to generalize efficient numerical integration formulas of the classical Chebyshev polynomials. 
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One of the most important application of the classical Chebyshev polynomials is the calculation of numerical 
quadratures - formulas equating a weighted integral of a polynomial function not exceeding a specific degree 
with a linear combination of polynomial values at some points called nodes 5p5 . The main point of quadrature 
formulas, which are in a multidimensional setting known as cubature formulas, is to replace integration by 
finite summing. The specific degree, to which cubature formulas hold exactly for polynomials, represents a 
degree of precision of a given formula. 

If a cubature formula holds only for polynomials of degree at most 2d — 1, then the least number of nodes 
is equal to the dimension of the space of multivariate polynomials of degree at most d — 1. Any cubature 
formulas which satisfy the lowest bound of the number of nodes have the maximal degree of precision and are 
called the Gaussian cubature formulas [4j[l7|[T8j[26 . In addition, it is known that such formulas exist only 
if the number of real distinct common zeros of the corresponding orthogonal polynomials of degree d is also 
exactly the dimension of the space of the multivariate polynomials of degree at most d — 1 ( 4 ] . It appears 
that among the 16 types of cubature formulas, resulting from 16 discrete transforms, 4 are Gaussian. The 
remaining 12, even though they are not optimal, extend the options for numerical calculation of multivariate 
integrals. 

Both antisymmetric and symmetric cosine functions can be restricted to their fundamental domain - a 
certain simplex in M n . For practical applications, it is essential that action of the permutation group S n on 
this simplex results in the entire space M n . Thus, a domain of any function of n real variables can be split 
into blocks - covered by sufficiently small copies of the fundamental domain. Since the fundamental domain is 
a subset of the n —dimensional cube, this method results in n\ times more blocks than the standard splitting 
into n —dimensional cubes. The multidimensional (anti)symmetric DCTs and the direct calculation of the 
corresponding interpolations are then performed in each of the blocks separately. Moreover, the resulting DCT 
coefficients can be further analyzed and serve as an input to multidimensional analogues of data hiding 19 
and image recognition 21 methods. Even though these methods can be applied to any sampled function, 


the more a given data set obeys the (anti)symmetries and boundary conditions of the DCTs at hand, the 
more suitable is the given method. The applicability of these methods to interpolation of potentials 


29 


can 

be tested - for instance, any potential in physics, which depends on the radial distance from the origin only 
(gravitational, Coloumb, Yukawa), is also symmetric as a function in M 3 . 

The 16 developed cubature formulas offer new options for numerical integration on the transformed funda¬ 
mental domain via a certain substitution. The transformed domain is, however, of a nonstandard shape - two 
main approaches, discussed together with applications in discretizations of PDEs and spectral approximations 
in [23 ,[24], can be followed to rectify this impracticality. Both these methods result in integration formulas on 
the simplex of the shape of the original fundamental domain. The first method straightens the transformed 
domain back to the original via a suitable straightening map; this methods results in higher densities of the 
nodes around the corners and edges of the simplex. The second method inscribes the original fundamental 
domain into the transformed domain and considers only the intersection points as nodes. The main advan¬ 
tage of both methods is that the resulting discrete calculus on the simplex allows a block decomposition of 
any bounded domain and thus, similarly to the (anti)symmetric cosine functions interpolation methods, can 
process practically any data set. 

In Section [ 2 J notation and terminology as well as relevant facts from 11 are reviewed. Simplified forms 
of three special cases of the generalized cosine functions, which appear as denominators in the definition 
of the Chebyshev-like multivariate polynomials, are deduced. The continuous orthogonality relations of the 
generalized cosine functions are presented. In Section [3J the one-dimensional DCTs of types V-VIII are 
reviewed and generalized to the antisymmetric and symmetric multivariate DCTs. The interpolation formulas 
in terms of antisymmetric and symmetric cosine functions are developed. In Section [4j the four families of 
the multivariate Chebyshev-like polynomials are introduced and three-dimensional examples presented. It is 
proved that each family of polynomials forms an orthogonal basis of the space of all polynomials with respect 
to the scalar product given by a weighted integral. In Section [5| the cubature formulas for each family of 
polynomials are derived. The last section contains concluding remarks and addresses follow-up questions. 


2. Symmetric and antisymmetric multivariate cosine functions 


2.1. Definitions, properties, and special values. 

The symmetric and antisymmetric multivariate generalizations of the cosine functions are defined and their 
properties detailed in 11 . The antisymmetric cosine functions cos^(x) and the symmetric cosine functions 
cos^(x) of variable x = (xi,...,x n ) G M n and labeled by parameter A = (Ai,...,A n ) G M n are defined 
as determinants and permanents, respectively, of the matrices with the entries cos(7rA iXj), i.e., taking a 
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FIGURE 1. The contour plots of the graph cuts (x% = |) of the antisymmetric cosine functions 
cos^(x). The cut of the boundary of the fundamental domain is depicted as the black 

triangle with the dashed lines representing the part of the boundary for which all the antisymmetric 
cosine functions vanish. 


permutation a G S n with its sign sgn(cr) one has 

COS^ (x) = T] S S n U) COS U-Vl^l) COS (7rA (T ( 2 )X 2 ) ■■■COS (tt\ a ( n )X n ) , 


CreSn 


cos+(x) = E cos {'k\ (T {i)Xi) cos (7 rA cr ( 2 )X2) • • • cos (7rA a ( n )X n ). 
cres n 


( 1 ) 


( 2 ) 

( 3 ) 


The cos A (x) functions are continuous and have continuous derivatives of all degrees in M n . Moreover, they 
are invariant or anti-invariant under all permutations in a G S n , i.e., 

cos“00)) = Sgn(cr) cos”0), cos“ (A) 0) = sgn(cr) cos” (x ), 

cost (°V)) = C0S A ( X ) ) C0S t(A) ( X ) = C0S A ( x ) > 

where a(x) = x ct(2 ), ..., a: CT („)) and a( A) = (A ct( i), A ff ( 2 ),..., A CT ( n) ). Additionally, they are symmetric 

with respect to the alternations of signs - denoting r* the change of sign of Xi, i.e., ri(x) = r^Oi, * • •, x u ..., x n ) = 
Oij • • • 5 ~ Xi ,..., x n ) it holds that 

cos“ (Ti(x)) = cos” (x) , cos“ (A) (x) = cos“ (x ), 

cost ( T i ( X )) = C 0 S f( X ) > cos t(A) ( a: ) = C0S A ( X ) ■ 

Considering the functions cost (x) with their parameter having only integer values k £ Z" and denoting 


( 4 ) 


P = 


1 1 

2 ’ ’ 2 


( 5 ) 


the symmetry related to the periodicity of the one-dimensional cosine functions, i.e., for t — (ti,... ,£ n ) G Z' 
we obtain 

CQSfc (x + 2 1) = cos*; (x), cos^ +p (x + 2 1) = cos*; +i0 (x), 


COS 

\t lH-h t n 


COS 




(A 


(6) 


cos+ ( x + 2t ) = cos fe 0). cos+ +p (x + 2i ) = (-!) 

Introducing the two index sets 

P+ ={(k 1 ,...,k n )GZ n | fci >fc 2 > ••• >k n > 0}, (7) 

P++ = {(fci,..., k n ) £ Z n I fc! > fc 2 > • • • > k n > 0} , (8) 

the relations {[2]), ((3j|, and @ imply that we consider only the following restricted values of k G Z n , 

cos*: (x), cos*T +p (x) : k £ P ++ , 

COS 


-+(x),cos+ (x) : keP + . 


( 9 ) 


By the relations we consider the antisymmetric and symmetric cosine functions cos^(x) labeled by 

k G Z n on the closure of the fundamental domain ) of the form 

F(S^ ff ) = {(xi,x 2 , ..., x n ) G M n I 1 > Xi > x 2 > • • • > X n > 0} . (10) 

Due to <®> <§, and the identity cost (ki + |) =0, valid for ki G Z and i G {1,... , n}, we omit those 
boundaries of F(S^ ff ) for which 

• — ^2+1, i G {!)••• , n — 1} in the case of cos^, (x), 


• Xi = Xi+ 1, i G {1,..., n — 1} or x\ = 1 in the case of cos fe+ (#), 

• = 1, i G {1,..., n} in the case of cos^ +p (x). 
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FIGURE 2. The contour plots of the graph cuts (#3 = |) of the symmetric cosine functions cos^ (x). 

The cut of the boundary of the fundamental domain F(5'f ff ) is depicted as the black triangle. 

The contour plots of the cuts of the three-dimensional antisymmetric and symmetric cosine functions are 
depicted in Fig. [l] and Fig. [2] respectively. 

For several special choices of k G Z n , the functions ([l]) are expressed as products of the one-dimensional 
cosine and sine functions | lT] . In addition to the formulas from ll], we calculate this form for the special 
cases needed for the analysis of the related orthogonal polynomials. Denoting 

pi ={n- l,rc — 2, 


3 1 


P 2 = \n- -,n- -, 


we derive the form of the functions cos+(x), cos pi (x) and cos p2 (x) in the following proposition. 

Proposition 2.1. Let k G N is given by 


^ 2 ^ for n odd, 
§ for n even. 


Then it holds that 

cos-(zi,...,x„) = (-l)^- 1 ) 2 0 s in (|(x* + x j)j sin (|(x; - Xj)j 




P 2 


cos; 


n 

(xi,... ,x n ) — (-l) fe 2 n(n_1) 0cos 0 sin (^(x* + xj)J sin (^(x* - x.,-)) , 

i=1 ' " ‘ ' 

n 

nip cos (|xj) • 


l<z<j<n 


i= 1 


where pi, p 2 , and p are given by ( [T2| ), and 0, respectively. 
Proof. We have from definition 0 that 


cos (xi,... ,x n ) = det 


/ cos (7r(n—l)cci) cos (7r(n—2)cci) ••• cos(7rcci) 1" 
cos (7r(n— 1 )^ 2 ) cos (7r(n— 2 ) 2 : 2 ) cos ( 7 ^ 2 ) 1 


\ cos (7r(n—l)x n ) cos (7r(n—2)a: n ) ••• cos(7rcc n ) 1 / 

Using the trigonometric identity for the powers of cosine function 

2 m-l CQS m q _ f cos {mO) + (T) C° S ( m - 2k ) 8 if m is odd > 

fcos ( [md ) + \ (m) + Efc=i 1 (U cos ( m ~ 2 k)9 if m is even, 
we bring the determinant to the following form: 


cos pi (xi,...,x n ) = (-l) fc 2 2 


0 det 


/I COs(7TfCi) ••• COS n 2 (7TfCi) COS n 1 (iTXi) ' 
1 COS ( 7 TX 2 ) ••• COS n_2 (7TCC2) COS™ -1 ( 7 ^ 2 ) 


( 11 ) 

( 12 ) 


(13) 

(14) 

(15) 

(16) 


\ 1 COS (7 TX n ) ••• COS™ 2 (7TCC n ) COS n * (7 TX n ) , 

Taking into account that the last determinant is of the Vandermonde type and that the trigonometric identity 

cos ( 7 TXj) — cos ( 7 TXi) = 2 sin + x j)^J sin (~^( Xi ~ x j )) 


holds, we obtain that 


(n-2) (rt -l) 


cos (xi,... ,x n ) = (—1)^2 2 


| | (cos ( 7 TXj) — COS (jTXi)) , 


l<i<j <n 


and the identity (14) follows. The proof of formula (15) is similar and relation (16) follows directly from 
definition §. □ 
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Note that the equalities (14)—(15) allow us to analyze the zeros of the functions cos pi , cos p2 , and cos+: 

• cos” {x) = 0 if and only if x = (aq,..., x n ) G F{S^) satisfies Xi = Xi+ i, 

• cos p2 (x) = 0 if and only if x — (aq,..., x n ) G F{S ^ ff ) satisfies Xi = aq+i or aq = 1, and 

• cos+(x) = 0 if and only if x = (aq,..., x n ) G satisfies aq = 1. 

Observing that all zero values are located on the boundaries of F{S ^ ff ) we conclude as follows. 

L4.J.J J-iiv, lUllUUlUllO 5 + 

F(5f) 


Corollary 2.2. The functions cos pi , cos p2 , and cos+ are nonzero in the interior of the fundamental domain 


2 . 2 . Continuous orthogonality. 

The antisymmetric and symmetric cosine functions are mutually continously orthogonal within each 
family. Denoting the order of the stabilizer of the point k G M n under the action of the permutation group S n 
by H k , i.e., 


H k = # {ak = k | a e S n }, 

and introducing the symbol h k = h kl • • • h kri by 


(17) 

h _ f1 if ki = 0, 

1 \ otherwise, 



we obtain 



/ cos^ (x) cos^, (x) dx = h k 5 kk > , 

Jf(s^) 

k, k! e P ++ , 

(18) 

/ _ GOSfe + (a;) cos*, + (a;) dx = 2~ n 5 kW , 

fc, k! € P ++ , 

(19) 

/ cos£ (x) cos£, (x) dx = h k H k 5 kk ’, 

k, k’ e P+, 

(20) 

/ cosL (x) cosY (x) dx = 2~ n H k S kk ’ , 

iF(Saf f ) 

e P+. 

(21) 


The orthogonality relations (18), (20) are deduced in 11 from the continuous orthogonality of the ordinary 
cosine functions cos (irmO ), 


/ cos (nmO) cos {irm! 0) = h rn 5 rnrn ^ m, m G Z . 

Jo 


The remaining two orthogonality relations follow from the continuous orthogonality of the cosine functions 
cos ( 7 r (m + 1/2) 0 ), 


/■ 


cos ( 7r ( m + - ) 0 ) cos ( 7r [ m' + - ) 0 ) = -h, 


1 


1 


( 22 ) 


Since the antisymmetric cosine functions are anti-invariant with respect to all elements of S n and the group 
S n applied on F(S^f) gives the cube [0, l] n , we have 

If gaff C0S Up( x ) cos k'+p ( x ) dx = I J oin cos k + P ( x ) cos k'+ P ( x ) dx • 


Using the property eos /c+p (x) = cos a ,^ k+p Ja'(x)) for a' G S n , which follows from Q, together with definition 
([!]), we have 


cos. 


°k+p( X ) COS k>+p( X ) = J2 Sgn U' 7 ') II COS ^ 3V(i)) COS ^7T (fc' + ^ 

Therefore, we rewrite the integral © 


V<j' (z) 


as 


J^_ ff) cos fc+p (x) cos fe , +p (x) dx = ^ sgn (cr) / cos ^tt (fc CT(i) + 0 x^ cos ^tt (fc' + 0 x^) dx t , 


where we made the change of variables from av(i) to aq. Using the relation ( [22] ) , we obtain 


X 




n 1 

cos“ +p (x)cos“ +p (x)dx = ^ sgn (cr) JJ -6, 


a es n 


i= 1 


2 ” ^cr(i) ^ 


Since k a (^ = for alH G {1,..., n} if and only if cr is the identity permutation and hi = k[, the orthogonality 
relation (19) is derived. The orthogonality relation (21) is obtained similarly. 
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3. Discrete multivariate cosine transforms 


3.1. DCTs of types V-VIII. 

The one-dimensional DCTs and their Cartesian product generalizations have many applications in math¬ 
ematics, physics, and engineering. They are of efficient use in various domains of information processing, 
e.g., image, video, or audio processing. All types of DCTs arise from discretized solutions of the harmonic 
oscillator equation with different choices of boundary conditions applied on its two-point boundary [3,27 


These conditions moreover can be applied at grid points or at mid-grid points which are distributed uniformly 
in the one-dimensional interval. Requiring the Neumann conditions at both points of the boundary at grid 
points generates the transform DCT I; the Neumann conditions at both points applied at mid-grid generate 
DCT II. Requiring the Neumann condition at one point and the Dirichlet condition at the other point at 
grid points generates the transform DCT III; the Neumann condition and the Dirichlet conditions applied 
at mid-grid generate DCT IV. Note that for these four types of DCTs both boundary conditions are always 
applied simultaneously at grid points or at mid-grid - the other DCTs of types V-VIII are obtained similarly 
by applying one condition at grid point and the other at mid-grid. The DCTs express a function defined on a 
finite grid of points as a sum of cosine functions with different frequencies and amplitudes. They have numer¬ 
ous useful properties as the unitary property and convolution properties [3]. The multivariate symmetric and 
antisymmetric generalizations of the cosine transforms of types I-IV are found in [9, 11 . In order to derive 
the symmetric and antisymmetric DCTs of types V-VIII, which have not so far appeared in the literature, 
we first review their one-dimensional versions. 


3.1.1. DCT V. 

For TV E N are the cosine functions cos ( irks ), k = 0,..., TV — 1, defined on the finite grid of points 


5 e 


2r 


27V — 1 


I r = 0,..., TV — 1 


pairwise discretely orthogonal, 

N—l 


E 

r=0 


2irkr 
27V — 1 


2nk'r 
2N - 1 


2N — 1 

4c& 




where 


\ if r = 0 or r = V, 
1 otherwise. 


(23) 


(24) 


(25) 


Therefore, any discrete function / given on the finite grid (23) is expressed in terms of cosine functions as 

N—l A N—l 


f(s ) = ^2 A k cos (irks), A k = 


k =o 


4 c k 

27V-1 


E c -/ 


r=0 


2 r \ ( 2 tt kr 

cos 


2N — 1 


2N — 1 


(26) 


The formulas (26) determine the transform DCT V. 


3.1.2. DCT VI. 

For TV E N are the cosine functions cos (irks), k = 0., N — 1, defined on the finite grid of points 

f 2 (r + 4) . 1 

s € < V V r = 0, ...,7V — 1> 


27V- 1 


pairwise orthogonal, 


N—l 

E 

r=0 


C r+ 1 COS 


2irk (r + 
2N- 1 


'2t rk' (r + ±) \ _2N-1 


2N-1 


4 Ck 




(27) 


(28) 


where c r are determined by (|25j). Therefore, any discrete function / given on the finite grid (|27j) is expressed 

cos | ——^I • (29) 


in terms of cosine functions as 
N—l 


/( s ) = E Ak cos ( 7rfcs )’ Ak = 2 tv"— 1 E Cr + 1 f 


N—l 


'2 (*■+§)' 


fc =0 


r=0 


27V — 1 


27V — 1 


The formulas (29) determine the transform DCT VI. 
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3.1.3. DCT VII. 

For TV G N are the cosine functions cos (7r (fc + ^) s), k = 0,..., TV — 1, defined on the finite grid of points 
(23), pairwise discretely orthogonal, 

N—l 


c r cos 


r=0 


27r (fc + |) r" 
27V- 1 


cos 


27T (fc' + \) r\ 27V - 1 
27V-1 ) = 4c fc+ i 


dkk' , 


(30) 


where c r is determined by (|25j). Therefore, any discrete function / given on the finite grid (|23j) is expressed 

N—l 


in terms of cosine functions as 

N—l 


f(s) = 57 Acos (tt (fc+ |) s), A 


k=0 




27V- 1 


r=0 


2 r 


27V — 1 


cos 


27t (fc + 


27V- 1 


(31) 


The formulas (31) determine the transform DCT VII. 


3.1.4. DCT VIII. 

For TV G N are the cosine functions cos (it (k + |) 5 ), k = 0,..., TV — 1, defined on the finite grid of points 


5 G 


2 ( r +1) 


2N + 1 


r = 0. N—l 


(32) 


pairwise orthogonal, 


N—l 

E 

r=0 


COS 


2A^+ 1 


COS 


' 2tt (fc' + l) (r+|) \ = 2N- + 1 
2A^ +1 I 4 




(33) 


Therefore, any discrete function / given on the finite grid ([32]) is expressed in terms of cosine functions as 


N-l 


k=0 


/( s ) = E A cos (7T (k+ |) s), A k = 


N—l 


27V + 1 


E/ 

r=0 


2N + 1 


COS 


2?r (k+ |) (r + |) 


2AT + 1 


(34) 


The formulas (34) determine the transform DCT VIII. 


3.2. Antisymmetric discrete multivariate cosine transforms. 

The antisymmetric discrete multivariate cosine transforms (AMDCTs), which can be viewed as antisym¬ 


metric multivariate generalizations of DCTs, are derived using the one-dimensional DCTs from Section 3.1 


The four types of AMDCT, connected to DCTs I-IV, are contained in [11]. Our goal is to complete the list of 
AMDCTs by developing the remaining four transforms of types V-VIII. First, we introduce the set of labels 


Dn = {(ki,... ,k n ) G Z n | TV - 1 > ki > 0, i = 1,... ,n} 
and to any point (fci,..., k n ) G Djv we assign two values 

dk = Ck x • • • c kn , 
dk = c/d+i • • • Cfc n +i, 


(35) 

(36) 

(37) 


where c/ Ci are determined by (25). The discrete calculus is performed on three types of grids inside 
two of them are subsets of two cubic grids, 


Cn = 
Cn = 


2n 


2 r n 


27V — 1 ’ ’ 27V — 1 

_ f /2 (7*1 + 4 ) 2 ( r n + \) 


27V — 1 


27V — 1 


Or, ...,r n )eD N 

(ri, ...,r n )eD N \. 


(38) 

(39) 


To any point sGdjv, which is according to (38) labeled by the point ( 7 * 1 ,..., r n ) G Dn, we assign the value 

=C ri . ..Cr n , (40) 


and to any point s G Cjv? which is according to (39) labeled by the point ( 7 * 1 ,... ,r n ) G ZVv, we assign the 
value 


S s — C ri _|_i . . . C rn _|_i. 


(41) 


In the following subsection we detail the proof of the AMDCT V transform. The proofs of the remaining 
transforms VI-VIII are similar. 
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3.2.1. AMDCT V. 

For TV G N we consider the antisymmetric cosine functions cos^T ( s ) labeled by the index set 


and restricted to the finite grid of points contained in Cat C 

F V- IV 2y h 2r n \ I / 

N - \\2N- i 7 '” 7 27v-1y 1 

The scalar product of any two functions /, g ■Fn~ -A M given on the points of the grid Fn~ is defined by 


> k n } 

(42) 

) e • 

(43) 


(. f,g)= 57 £ sf( s )9( s )’ 


(44) 




where e s is given by (40). Using the orthogonality relation of one-dimensional cosine functions (24), we 
show that the antisymmetric cosine functions labeled by parameters in &, k' G are pairwise discretely 
orthogonal, i.e., 


(cos k , COS fe , ) = 57 C0S fe COS k> ( S ) = d k 1 ( 27V 4 M 


where dk is given by (36). Denoting the set of points 

= eCw 


ri > r 2 > • • • > r. 






(45) 


(46) 


we observe that the functions cos^, ( s ), k G vanish on the points seF N \ F%’- and thus 

(cos“,cos“)= 57 £ s h 7 1 cos k ( S ) COSfe, (s), 


sEFjy 


with the order of the stabilizer i7 s determined by Acting by all permutations of S n on the grid Fn we 

obtain the entire cube Cat, i.e., S n Fjy = C/v, and therefore we have 


N—l 


(cos k ,cos k ,) = 4 ^ 


C r 1 • • • c r n COS/c 


ri,...,r n =0 


2ri 


2r„ 


27V — 1 7 27V - 1 


cos,, 


2n 


2r n 


27V — 1 7 27V - 1 


Moreover, the product of two antisymmetric cosine functions labeled by &, k’ G D N is rewritten due to the 
anti-invariance under permutations as 


cos^, (x) cos k , ( x ) = ^ sgn (aa f ) JJ cos ( 7 r& aa /(q£ a /(q) cos ( 7 r^-V(i)) 
a,a'ES n i= 1 

Together with the S'n-invariance of the set ZAv this implies that 

n N—l 


, _ _v i r / 27T k a (i}Ti\ ( 27rk f -Vi 

{cos k ,cos fc ,) = 2_^ sgncr j j 2^ Cr * cos * I cos 

creS n i= 1 n=0 


27V — 1 


27V — 1 


Finally, we apply the orthogonality relation of one-dimensional cosine functions (24) to obtain (45), 

j-i /2-/V — 1' ~ 


_ 27V — 1 

(cos k , COS^,) = 57 s S n ( <T ) n 17- S **<.*K = ^ 

<T gs„ i =1 




Due to the relation (45), we expand any function / : F^’ 


in terms of antisymmetric cosine functions 


as follows: 


/( s ) = 57 A COS fc (s) with A = ( 2iV 4 _ 1 ) 57 e */( s ) cos fc ( s ) • 


keD~ 


?n’~ 


Validity of the expansion follows from the fact that the number of points in D N is equal to the number of 
points in F^~ and thus the functions cos^ with k G D ^ form an orthogonal basis of the space of all functions 


S-Fi 


. with the scalar product (44). The remaining transforms AMDCT VI-VIII are deduced similarly. 
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3.2.2. AMDCT VI. 

For iVeNwe consider the antisymmetric cosine functions cosjjT(s) labeled by the index set k G D ^ and 
restricted to the finite grid of points contained in Cn C F(S ^), 


^ VI - = 
In — 


( 2 ( r i +1) 2 ( r » +1) ) 

y 2N-1 2N-1 J 


(ri,.-.,r n ) 



(47) 


The antisymmetric cosine functions labelled by parameters in k, k' G D N , are pairwise discretely orthogonal 
on the grid F^ 1, ~, i.e., 


E e a cos k (s)cos k ,(s) = d k 1 

^ f n'~ 


/ 2N — 1 
V 4 


n 


(48) 


where e s is given by (41). Therefore, we expand any function / : Lv 
functions as 


VI,- 


M in terms of antisymmetric cosine 


f(s) 


( 4 \ 

E A k cos k (s) with A k = d k [ 2jV _ - J Y £ V(s) cos k (s). 

keD~ s ^ f n 1 ’~ 


3.2.3. AMDCT VII. 

For TV G N are the antisymmetric cosine functions cos^^s) , k G D and p = (4,... 5 4) restricted to the 

finite grid of points = F^~ C Cn pairwise discretely orthogonal with respect to the scalar product 

(44b, i.e., 


E £ * cos k + P ( s ) cos fe'+ P ( s ) = d k 1 ( 2iV 4 - ) 

S e^ n ’- V 7 


where dk is given by (37). Therefore, we expand any function / ■■Fn 
functions as follows: 


VII,- 


$kk' i (49) 

in terms of antisymmetric cosine 


m = E A - COS k+p (s) with A k = d k l^n) E e s /(s) cos fe+p (s) 

V 7 seiV 1 ’ - 


3.2.4. AMDCT VIII. 

For TV G N are the antisymmetric cosine functions cos (s) , k G and p = (4, ..., |) restricted to the 
finite grid of points 


r^VIII- 

■ AT 


2 ( r l + |) 2 ( r n + !) 


2JV + 1 


2iV+ 1 


(n,...,r„) e D 


N 


(50) 


pairwise discretely orthogonal, i.e., for any &, k' G it holds that 

Y cos k + P ( s ) cos fe'+ P ( s ) = ( 2iV . + 1 ) <w • 




VIII,- 

N 


(51) 


Therefore, we expand any function / : -A M in terms of antisymmetric cosine functions as follows: 

/(-)- E Acos- +p(s ) with A k =( 4 ) JZ /(s) cos- +p (s). 

1 7 7 \ Z ' 1 / _ 7-,VIII.- 


keD~ 


seF, 


VIII,- 

N 


3.3. Interpolations by antisymmetric multivariate cosine functions 


Suppose we have a real-valued function / given on F(S^). In Section 3.2 we defined three finite grids in 


F(S^), namely, F^ ’ , = F^ ’ , and F^ ’ . We are interested in finding the interpolating polynomial 

of / in the form of a finite sum of antisymmetric multivariate cosine functions labeled by k or by k + p, k G D ^ 
in such a way that it coincides with / on one of the grids in F(S^). We distinguish between four types of 
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Figure 3. The cut of the model function (53) with fixed parameters (54) and z = 




Figure 4. The antisymmetric cosine interpolating polynomials i/’w |) of the model 

function (53) shown in Fig. [3] with N = 5,10,15. 

interpolating polynomials defined for x G and satisfying different conditions: 


'n (• X ) = E B kCOS k {x), 

(«) = /(«), sgF 

keD~ 


’n~( x ) = E B kCOS k+p (x), 

= sgf y-, 

ke-D- 


W n ’~( x )= V B kCOS~(x), 

= /( S ) > S G F N h ~ . 

keDJ, 


-r •-(*)= v B k cos k +P ( x )> 

Vv ( 5 ) = f( s ) > 5 ^ 


(52) 


keD: 


According to Section [372] the coefficients Bk are chosen to be equal to corresponding A^. In fact, it is not 
possible to have other values of the coefficients Bk since it would contradict the fact that the antisymmetric 
cosine functions labeled by k or by k + p with k G D^ are basis vectors of the space of functions given on the 
corresponding grid. 

Example 3.1. For n = 3 we choose a model function / and interpolate it by (x,y, z) and y, z) 

for various values of N. The model function / is a smooth function given by 

(x - x 0 ) 2 + (y- y 0 ) 2 + {z- Z 0 ) 2 ) 


fa,p,(xo,yo,zo)(x,y,z) = exp 


2a 2 


+ P 


with fixed values of the parameters <a, /3, (xo, yo, zq) as 

= 0.079 , p = 3 , (x Q , y 0 , z 0 ) = (0.8,0.54,0.3) 

Figure [i] 


(53) 

(54) 


contains the function / for fixed parameters (54) in the fundamental domain F(S'| ff ) and cut by the 
plane z = \. The antisymmetric interpolating polynomials of type V and VII for N = 5,10,15 are depicted 
in Figures |4] and [5j 

Integral error estimates of the approximations of the model function (53) by the antisymmetric interpolating 
polynomials of types V and VII, 


/„ 


f(x, y, z) - ipx’ ( x , y, z) 


are shown in Table [l] for several values of N. 


dx dy dz , 


L 


F(S * ff ) 


f{x,y,z) - ip™’ ( x ,y,z) 


dx dy dz , 













MULTIVARIATE COSINE TRANSFORMS AND POLYNOMIALS 


11 



3.4. Symmetric discrete multivariate cosine transforms. 

The symmetric discrete multivariate cosine transforms (SMDCTs), which can be viewed as symmetric 
multivariate generalizations of DCTs, are derived using the one-dimensional DCTs from Section 3.1 The 
four types of SMDCT, connected to DCTs I-IV, are contained in [ll]. Our goal is to complete the list of 
SMDCTs by developing the remaining four transforms of types V-VIII. Note that the coefficients e si s s , and 
dk/dk are given by ( |40| ), © and ( [36] ) 7 ( [37]) , respectively. In the following subsection we detail the proof of 
the SMDCT V transform. The proofs of the remaining transforms VI-VIII are similar. 

3.4.1. SMDCT V. 

For TV £ N we consider the symmetric cosine functions cos^ (s) labeled by the index set 

d n — {(^ 1 , ■ ■ ■ ,K) £ D N \ ki > k 2 > ■ ■ ■ > k n } ( 55 ) 

and restricted to the finite grid of points contained in Cjv C F(,S'.;) tr ), 

2 r n 


77>V,+ _ 

t n — 


\ \2N —I’’"’ 


2JV-1 


(n,. 


£ -Djv 


} 


(56) 


The scalar product of any two functions /, g : F^’ + —> R. given on the points of the grid F%' + is defined by 

(f,g)= e s Hr 1 f(s)g(s), f,g:F^’+^R, (57) 




where H s is given by ©• Using the discrete orthogonality of one-dimensional cosine functions ( [24] ) we show 
that the symmetric cosine functions labeled by parameters in k, k' £ D ^ are pairwise discretely orthogonal, 


i.e., 


cos^ , cos 


^2 £sH * lcos k( s ) c ° 4 '( s ) = 


H k ( 2N — 1 


seT 


.v,+ 

N 


$kk' 


(58) 


Note that acting by all permutations in S n on the grid we obtain the whole finite grid Cjv, where the 
points s £ Ftf, which are invariant with respect to some nontrivial permutation in S ni emerge exactly H s 
times. Since the symmetric cosine functions are symmetric with respect to *S n , we obtain 

N—l 


(cos+ COS+) = Ij 5Z 


+ ( 2ri 2 r n 

n\ ^ ° ri ’ " Crn C ° Sfe \2N - 1 ’ ’' ’ ’ 2N - 1 

ri,...,r n =0 x 


+ ( 2r i 

COSft ' 1 2N - 


2r„ 


2N-1 


where n\ represents the order of the group S n . Moreover, the product of two symmetric cosine functions 
labeled by k, k! G D~^ is rewritten due to the invariance under permutations as 

n 

cos^(x) cos)!) (x) = JJcos cos ( 7 rk-x^^) . 

cr,a'ES n i= 1 

Together with the ^-invariance of the set Dn this implies that 

/ + v - ^ T — r / 2'Kk a (i)Vi \ ( k'-Vi 

(cos£,cos£) = 2^ 11 cos 1 ^“ 1 cos 

a es n i=in=o 


2N-1 


2N - 1 
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Finally, we apply the one-dimensional orthogonality relation (24) to obtain (58), 


+ +\ \-TT 2N - 1 x H k f2N- 1 

<jes n i=i 


(cos,;, cos 






V 4 

Due to the relation ( [58] ) , we expand any function / : F^ ,+ -A M in terms of symmetric cosine functions as 

/(«) = A COS + (s 


/cGD 


A k = ^( 7^rz l ) H V00 cos fe 0 s 


if a- \2iV 


seF' 


Validity of the expansion follows from the fact that the number of points in D ^ is equal to the number of 
points in Fn’ + and thus the functions cos^ with k G D^ form an orthogonal basis of the space of all functions 
/ : F^’ + -A M with the scalar product (57). The remaining transforms SMDCT VI-VIII are deduced similarly. 

3.4.2. SMDCT VI. 

For TV G N we consider the antisymmetric cosine functions cos^ (s) labeled by the index set k G D ^ and 
restricted to the finite grid of points contained in Cn C 


^vi,+ 

** N 


2 (n + |) 2 (r„ + |) 


(ri,...,r„) e £$> • 


(59) 


The antisymmetric cosine functions labeled by parameters in k, k f G are pairwise discretely orthogonal 
on the grid F^ I,+ , i.e., 


Z) £ sH s 1 COS+ (s) COS+ (s) = —j— 


H k ( 2N — 1 




vi, + 

N 


k 


$kk' • 


(60) 


Therefore, we expand any function / : F^ I,+ -A M in terms of symmetric cosine functions as follows: 

/0) = Ak COS k( S ) with A k = jf ( 2 N- l ) ^ ^sH- 1 f(s) cos~ (s). 


keDt 


s 


vi, + 

N 


3.4.3. SMDCT VII. 

For V G N are the symmetric cosine functions cos^ +p (,s) , k G and p = ( 4 ,..., 4 ) restricted to the 

finite grid of points F^ II,+ = F^ ,+ C Cn pairwise discretely orthogonal with respect to the scalar product 
(57b, i.e., 


Z) £ sH a 1 COS + +p (s) COS+ +p (s) = y 


/ 27V — 1 




(61) 


Therefore, we expand any function / : F^ II,+ -A M in terms of symmetric cosine functions as 


N—l 


/(*)= 5Z A feCOS+ +p (s), -4fc = T ( 4 _ ) s sHPf{8) COS+ +p {s) . 

k ' 7 7,VII,+ 






3.4.4. SMDCT VIII. 

For V G N are the symmetric cosine functions cos^ +p (,s) , k G and p = ( 4 ,..., 4 ) restricted to the 
finite grid of points 


? viii,+ _ ( ^2 (tt + 4 ) 2 (r n + 4 ) 


A 


(Fi,...,r n ) G D+ 


11 27V + 1 5 ‘ ‘ ‘ 5 27V + 1 

pairwise discretely orthogonal, i.e., for any k,k' G D^ it holds that 

/ 27Y _j_ \ \ n 

V cos fc+p( s ) cos fe'+p( s ) = ( —4—) <w • 

Therefore, we expand any function / : F^ III,+ -A M in terms of symmetric cosine functions as 
/(«) = V Acos+ +p (s), A=T ^ 4 


H k \ 2N + . 


l) V VA cos+ +p (s). 


s€FjO iii,+ 


(62) 


( 63 ) 
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3.5. Interpolations by symmetric discrete multivariate cosine fu nctio ns. 

Suppose we have a real-valued function / given on F(S^). In Section 3.4 we defined three finite grids in 
F(S^), namely, ,+ , F\ y’ + = ,+ , and F^ ,+ . We are interested in finding the interpolating polynomial 

of / in the form of a finite sum of symmetric multivariate cosine functions labeled by k or by k + p, k G D^ 
in such a way that it coincides with / on one of the grids in F(S^). We distinguish between four types of 
interpolating polynomials defined for x G F(S^) and satisfying different conditions: 


F’ + ( X ) = E BkC OS k(x), 

keD+ 

4>n' + ( x ) = E BkCOS k+ P ( x )> 


Fn + {s) = f{s), S€F%’+, 


K l,+ F = f(s), seFj 


VI,+ 
N > 


keD^ 

E 

keD+ 


F n,+ ( x )= E B k cos+{x), 


F N U ’ + (s) = f(s), seF% n ’ + , 


(64) 


FF + f= E Bkcos i+ P ( x ^ 

keD+ 


F m ’ + (s) = f(s), seF N 


.VIII,+ 


According to Section |3i4} the coefficients B & are chosen to be equal to corresponding A^. In fact, it is not 
possible to have other values of the coefficients B & since it would contradict the fact that the antisymmetric 
cosine functions labeled by k or by k + p with k G D ^ are basis vectors of the space of functions given on the 
corresponding grid. 


Example 3.2. For n m 3 is the smooth function (53) with values of parameters given by (54) chosen as a 
model function. We interpolate this function by polynomials of symmetric cosine functions ip N ’ + (x, y, z) and 
i/> N ,+ (x,p,z) with N = 5,10,15. Plots of this symmetric cosine interpolating polynomials are depicted in 
Figures [6] and [7| 






Figure 6. The symmetric cosine interpolating polynomials Vhv + ( x ,^, §) °f the model func¬ 
tion (53) shown in Fig. [3] with N = 5,10,15. 


Integral error estimates of the approximations of the model function (53) by the symmetric interpolating 
polynomials of types V and VII, 

2 r _ 2 

dx dy dz , 


f. 

yT 

1 

+ 

yT 

Th 

2 f 

dx dy dz , / 

yT 

1 

+ 

yT 

4E 






are shown in Table [l] for N = 5,10,15, 20, 25, 30. 

4. Chebyshev-like multivariate orthogonal polynomials 


Recall that the vectors p, pi , and p 2 are defined by and ( fl2| ), respectively. Introducing the 

functions Xi, X 2 ,..., X n given by 


Xi cos (i ? o,...,0) ’ X 2 COS (1,1,0,...,0) ’ COS (1,1,1,0,..., 




— + 


0) 


, ... , X n = cos+ A A) , (65) 
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Figure 7. The symmetric cosine interpolating polynomials '0^ II,+ (x, y,\) of the model func¬ 
tion (53) shown in Fig. [3] with N = 5,10,15. 


N 

fp(S* ££ ) f V’iv’ 

r n /VII,- 2 

Jf(s^) J Fn 

J‘f(s^ ££ ) f Vw’ + 

r _p /VII,+ 2 

JF(S^) J VAT 

5 

0.648691 

1.396870 

0.725031 

1.502161 

10 

0.007940 

0.007599 

0.007191 

0.006471 

15 

0.001350 

0.001407 

0.000440 

0.000492 

20 

0.001034 

0.001058 

0.000171 

0.000195 

25 

0.000835 

0.000847 

0.000084 

0.000097 

30 

0.000698 

0.000705 

0.000047 

0.000054 


Table 1 . Integral error estimates of the polynomial approximations of the model function / 
(53b for N = 5,10,15, 20, 25, 30. 


we demonstrate that the following defining relations, valid for all points from the interior of x E 

F(sf )°, k e P + , 


V I k ’ + (X 1 (x),...,X n (x))=cos+(x), 


V I k ’~(X 1 (x),...,X n (x)) = 


COS 


k+p 




cos Pl (x) 


V I k II ’ + (X 1 (x),...,X n (x)) = 


a k+p 


(x) 


COSp (x) 


r' II ’~(x 1 (x),...,x n (x)) = 


k~hp 2 


(x) 


( 66 ) 


determine four classes of orthogonal polynomials ,7^^’ E M[Xi,...,X n ] of degree fci. 

Since from Corollary 2.2 it follows that the three functions in the denominators of (|66|) have nonzero values 


cos P2 (x) 

’ ni ’-y k n e 


inside F(S^ ff ), the functions (66) are well defined for any point in the interior of F(S^). Moreover, whenever 
one of the denominators is zero at some points of the boundary of F(S^) the corresponding nominator is 
zero as well. We introduce an ordering within each family of polynomials ( [66] ) . We say that a polynomial p\ 
depending on k = (Aq,..., k n ) E P + is greater than any polynomial P 2 depending on k' = (k [,..., k' n ) ^ k if 
for all i E {1,..., n} it holds that ki>k[\ equally, we state that P 2 is lower than pi. Note that the functions 
(66) can be viewed as generalizations of Chebyshev polynomials of the first and third kinds [5]. 


4.1. Recurrence relations. 

The construction of polynomials is based on the decomposition of products of symmetric and antisymmet¬ 
ric cosine functions. There are three types of products which decompose to a sum of either symmetric or 
antisymmetric cosine functions. Such a decomposition is obtained by using classical trigonometric identities 
and is completely described by 


COst 

( A lr 

••’An) 

■ COS tu. 

= iy 

-,Mn) 2 n ' ^ 

E 

COS (Ai+ai^ cr ( 1 ),...,An+anMcr(n)) ’ 




cr es n 

r a^=±l 

\ 






i 

cos 7 

(Ai,. 

..,A n ) 

• cos E. 

- — V 

-,Mn) 2 n 

E 

^^(Ai+ai/i^/!) ,...,An+anMcr(n)) 5 




O-eSn 

r ai=± i 

\ 





\ i=l,...,n 

/ 

cos 7 

(Aiv 

..,A n ) 

' cos (m,- 

= ly 

-,Mn) 2 n ' ^ 

E 

sgn (cr) cost, , , , , 

& V ' (Al+ai/Xcr/!) ,...,An+0'nMcr(n)) 




creS n 

r ai=zLl 

} 


( 67 ) 
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Using (67), we obtain the following pertinent recurrence relations. Let A = (Ai,..., A n ), and let h be a 
vector in M n with ith component equal to 1 and others to 0, then 


cos a = l!(n — 1)! COS U A - cos f-2/ 1 - 2 (' 


COS^_; i+; . +COS^_; i _,.) , 


COS, = 


C0S A-Zi-i 2 A COS A-2Z 1 -2Z 2 C0S A-2i 1 COS A-2i 2 


A 2!(n — 2)! 

n 

- ( cos i-h+k +cos x-2h-h-u + cosf_ l2 _ h +c°s±_ 2h _ l2+li ) 


i =3 


n 

- E (■ 


COS A -h-h+k+lj +C0S A -h-h-k-lj + COS \-h-l 2 +li-lj + COS \-h-l 2 -li+l 


)■ 


= -7 cos±_ h _ l2 _..._ ln X n - ]T cos±_ 2/i - 53 < 

m) 




— • • • — cos 


\-2h-2h -2 l n 


( 68 ) 


i= 1 


Using the relations (68) and the symmetry properties of cos^ each polynomial of (66) is expressed as a linear 
combination of lower polynomials and a product of some lower polynomial with some Xi. Therefore, all 


polynomials (66) are built recursively. 


Proposition 4.1. Let k G P + . The functions V k ^ and V] 


-.hi, ± 


in variables Xi,... ,X n . The number of P k :± or with k\ = d is equal to the number of monomials of 


iii,± 


are expressed as polynomials of degree k± 


degree d, i.e., 


d + n — 1 
d 


(69) 


Proof. At first, we proceed by induction on k\ to show that any P^ ,+ is expressed as a polynomial of degree 


ki in the variables (65). 


• If = 0, P 7 ’+ is trivially the constant polynomial n\ of degree 0. 

• If k\ = 1, the polynomials Xi, X 2 ,... , X n of degree 1 correspond by definition to the set of functions 


? (1 k n V Where 1 ^ k 2 > ■ ■ ■ > k n > 


0. 


• If k\ = 2, using relations (68) and the basic properties of symmetric cosine functions, we deduce 

that any P 7 ^ - ^ k ) constructed from the decomposition of the products X{Xj. Indeed, we start 

by obtaining the lowest polynomial P 7 ^ 0 ) f rom the product X\ which decomposes into the linear 
combination of P 7 ^ 0 ) an d polynomials of degree 1. We continue by the decomposition of the 

product X 1 X 2 to have the polynomial expression of P 7 ^ and so on. Finally, we obtain that each 
P 7 2 t 2 k ) is expressed as a polynomial of degree 2 in variables (65). 

• We next suppose that every V I ^ 1 k2 k ^ is a polynomial of degree d— 1 in variables (65), then again 


using relations (68) the polynomial P 7 ^ 2 k ^ is a linear combination of lower polynomials of degree 


at most d and P 7,+ ~ ~ Xi of degree d. 

{d—l,k 2 ,...,k n ) 

By induction, this results in the fact that any P^ ,+ is expressed as a polynomial of degree k\ in the variables 
(65). By similar arguments, we obtain the same statement for the functions V k ~ and Pf 77,± . Note that we 


start with the products of cos pi , cos+, cos p2 , respectively, with the variables Xi. 

To prove the second statement, we observe that the number of polynomials of type T k ± or P^ 77,± with 
k\ = d is equal to the number of elements in {(& 2 , & 3 , • •., k n ) G Z n_1 | d > k^ > k% > • • • > k n > 0} which is 
the same as 


i=2 


H 1 ^ (^n + In-1 + ’ • * + h, • • • + In-li ^n) | hi ^0 £ Z~°, l() + l{ — d ^ 

or equivalently 

# |v a, A >2 • • • AN ai > 0, “ d \ = 

The proof of the last equality is found in ^. 


i= l 


d + n — 1 
d 


□ 
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Example 4.1. In particular, the relations (68) together with the symmetry properties of symmetric cosine 
functions imply the following set of recurrence relations for P^ k2 ks y 

^( 0 , 0 , 0 ) = ^ ’ ^( 1 , 0 , 0 ) = ^1 5 ^( 1 , 1 , 0 ) = X 2 ’ ^( 1 , 1 , 1 ) = ^3 , 


(70) 


/ci > 2, /c2 = fe = 0 : 

- 1 > k 2 > ks = 0 : 

ki — 1 > &2 = > 0 : 

- 1 > k 2 > k 3 > 0 : 

ki - 1 = k 2 > ks = 0 : 

k\ - 1 = k 2 > ks > 0 : 

ki — 1 = k 2 = &3 > 0 : 

k\ = k 2 > 2, & 3 = 0 : 


&i — A ^2 > ^3 + 2 > 2 


P( M) 0)-P^nn^ 


>h+ 

(/ci-1,0,0)" 


(fei —2,0,0) 


4,4 

(fei-1,1,0) ’ 


p 7 ’+ 

' (fei,fe 2> 0) 

"P, 


A^A 1 nAl - A'Ao n\ - 2A / y_ 1 ^ - A /,+ 


(fcl— 1,/C2,0)^ 


(fei-2,fe 2 ,0) 


(/ci-l,/c 2 ,l) 


(/Cl-1,/C2 + 1,0) 


A 


4,4 

(/ci-l,/e 2 -l,0) ’ 

4,4 _ T)k-\- 

— 7 (kx-IMM) 1 
4 5 + 

(fei — 1 ,/c 2 ,fe 3 ) 


(ki,k 2 ,k 2 ) 

1 ,4 

( k!,k 2 ,k 3 ) 


*i-A 


AA„_ = A,A , Wi - ‘Pk+n~ AA , , ,, - A 


4,4 

(ki—2,k 2 ,k 2 ) 

4 ,4 

(fcl 2 ,/u 2 ,fc 3 ) 


— 2P 7,+ — 2 / P /,+ 

/fc. _i (/ei-l,/c 2 ,/e 2 -l) ’ 


(/ei — l,/c 2 + l,fe 2 ) 

■4,4 _ q-yl ,-\- 

(fei-l,fe 2 ,fe 3 + l) ' (/Cl- 1 ,/C 2 ,/C 3 - 1 ) 


-A 


4 ,4 

(/Ci-l,/C2 + l,/c 3 ) 


-A 


4,4 


'p4 # *4 _ W _ T>k+ 

7 (fei,fei-l,0) — 2 7 (/ci-l,/ci-2,0) 


-p J *+ 

7 (fei-l,fei-l,l) ’ 


T)k+ _ ^LT> 7 ’ + V' _ T)k + 

' U.. ^ ' (U-.— 1 _ 1 T^^l ' I 


(fci,fci —l,fc 3 ) — 2 7 (^1 — Ijfei — 1 ^ 3)^^ 1 7 (/ci-1,/ci-2,/c 3 ) 2 (k 1 — l,k 1 — l,k 3 -\-l) 


_ Ip 7 ’+ 

2 7 (fei-i,fei-i,fe 3 -i) ’ 

p 7 ’+ _ 1^4,+ y _'p6+ 

7 (fei,fei-l,fei-l) — 3 7 (/ei —1,/ci —1,/ci —l)^ 1 7 (fei-l,fei-l,fei-2) 


^4,4 


P, 


4,4 


= 2P 7 ’ + 


(fei,fei,0) — -(/ci-l,,i-l,0)^ ^^ > (/b 1 -l,/c 1 -2,0)^ 1 Z/ ^i-l,fei-l,l)^ 1 


A,+ 


xXi - 2 A 


/,+ 


,X| 


+ A 


-^5 + 

(fei-2,fei-2,0) 


+ 4A 


^5+- 

(fei-l,fei-l,0) 


+ 4A 


1 5+ 

(fei-l,fei-2,l) 


+ 2A 


h + 

(fei-l,fei-3,0) 


+ 2A 


I, + 

Oi—1,/ci —1,2) ’ 


^l,fcl,fc 3 ) Z/ (/Cl— l,fci — 1,/C3)^ 2 


^^ + - 1 2 ^ > (/ c ^- 1 ,/ ci - 2 ,/ c 3 )^ 1 V - l ,^ i - l ,^ 3 + l) Xl 


— 'P 7 ^ 4- P 7 ’ + 4 - 2 P 7,+ 

7 (/ci-p/ci-p/cs-l)^ 1 ^ 7 (fei-2,fei-2,fe 3 ) ^ ^ 7 (fei-l,fei-2,fe 3 41) 


+ 2P. 




1 5 + 

(/Cl—1,/Cl—2,/c 3 —1) 

/,+ 


+ 4A 


-^5 + 

(fel —l,fel —l,fe 3 ) 


+ 2 P, 


/, 4 




/,4 


(fei-l,fei-l,fe 3 -2) ’ 


fci - fc 2 - k 3 + 2 > 3 : 2 ) 


(fei-l,fei-3,fe 3 ) ^ 7 (fci — 1,/ci — 1 ,/c 3 +2) 

2 . 


2A'; ! ,,._oA 2 -2pA_, ,._o,_,Ai - 4/A, t ,_4i 


Oi —1,/ci —1,/ei—2)" 


(fei-l,fei-2,fei-2)^ 


3 7 (fci — l,fci — 1 ,/ci — 1 )^ 


-A 


4,+ 


(/ci — 1,/ci —1,/ci —3) + ^ > (fei-2,fci-2,fei-2) 


/,+ 


+ 5P ( 


h + 

(fci — l,fci — l,fci— 2) 


+ 4A 


/ci = fc 2 = ^3 + 1 > 2 : A 


/,+ 


-^5 + 

(fei — l,fci — 2,fei—3) 

- 


+p, 


-^5 + 

(k%— l,fci — l,fci— 4) ’ 


(fci ,/ci ,/ci —1) 3 7 '(/ci-l,/ci-l,/ei-l)^ 2 ^ > (/ci-l,/ci-l,/ci-2) y ^ 1 r (fei-l,fe 1 -2,fei-2) 


/,+ 




4 5+ 


+ A 


fci = &2 = &3 > 3 : 




(fei",fei,fei) — o 1 ,/ci — 1 ,/ci — l )^ 3 6 V_ 1 ,, 1 _ 1 ,, 1 — 2 )X 2 + 3 V_i,,i- 2 ,,i- 2 )X 1 


-^5 + 

(fci —l,fei —l,fei —1) 


+ A 


4,4* 

(fci-l,fci-l,fci-3) ’ 


4,4 


+ 2 P ( 
-9P ( 


/,4 


(fei — l,fei — l,fei —1) ^4 


Xi + 3)P ( 


4,4 


Xi — V 

(k 1 — l,k 1 — l,k 1 —3) yY l 7 (fei-2,fei-2,fci-2) 


4,4 


4,4 


6 P 7,+ 

U/ Oi —1,/ci—2,/ci—3) 


3P 7 ’+ 

07 (fei-l,fei-l,fei-4) 


(71) 


(/ci—1,/ci —1,/ci—2) 

together with additional polynomial expressions 

P ( J ’+ 0 ) =2X| - 2X( - ^XiX 3 + 8 X 2 + 6 , 

A 7 ri,i) = l X2X 3 - AX 2 + 2 X 3 + X,, 

?mr\ x l - 6X 2 + 3Xf + 4X a X 2 - 12X 2 - 6 , 

i)”|x|X 3 - «X!Xf - 2X 2 X 3 + 8 X 2 X 3 + 9 X 3 , 

44,3) = y ^3 - 12X|X 3 + 8 X 1 X 3 2 + 9X?X 3 - 36X 2 X 3 - 27X 3 . 

Similarly, one may find recurrence relations for 'P( k ~ k2 k y and P 7 ^ 7 ’^ k ^. The polynomials 'P( k f k2 k y and 
P 777 ,t feg ) of degree at most two are shown in Tables |2j—1^| 


(72) 
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1 

* 2 

*1 

* 1*3 

*2 

*3 

p J ’+ 

7 (0,0,0) 

6 






p 1 ’+ 

7 (1,1,0) 

0 

1 





T> 7 ’ + 

(2,0,0) 

-6 

-4 

1 




p J ’+ 

' (2,1,D 

0 

-1 

0 

1 

3 



T> 7 ’ + 

(2,2,0) 

6 

8 

-2 

4 

3 

2 


T> 7 ’ + 

(2,2,2) 

-6 

-12 

3 

4 

-6 

4 

3 



1 

* 2 

*T 

* 1*3 

*1 

*3 

•pi,- 

7 (0,0,0) 

1 






pi 

7 (1,1,0) 

3 

2 





T>u- 

(2,0,0) 

-4 

-2 

1 




pi,- 
7 (2,1,1) 

-3 

-2 

1 

2 

3 



(2,2,0) 

12 

14 

-3 

4 

3 

4 


(2,2,2) 

-9 

-12 

3 

10 

3 

-4 

8 

9 


Table 2. The coefficients of the polynomials P* k f k2 ks ) with < 2 and ki + k^ + k% even. 



V 

*3 

* 1*2 

* 2 * 3 

T > 7 ’- 

(1,0,0) 

1 




7 (TUI) 

1 

2 

3 



'pi , — 

7 (2,1,0) 

0 

4 

3 

2 


'D 1 ’- 

7 (2,2,1) 

1 

8 

3 

0 

4 

3 



V 

* 3 

*1* 2 

*2*3 

7 (i,o,o) 

1 




7 (EM) 

0 

1 



Tp ’ + 

7 (2,1,0) 

-1 

-1 

1 

2 


T> 7 ’ + 

7 (2,2,1) 

1 

2 

-1 

2 

3 


Table 3. The coefficients of the polynomials P( k ± k2 k3 ) with k\ < 2 and k\ + k^ + odd. 



1 

*i 

*2 

*3 

*i 

*1*2 

*1*3 

*2 

X 2 X 3 

*3 

pin,+ 

7 (0,0,0) 

1 










plll>+ 

7 (1,0,0) 

-1 

1 

3 









plll>+ 

7 (1,1,0) 

1 

2 

3 

2 

3 








p 111 ■> + 

7 (1,1,1) 

-1 

1 

-2 

4 

3 







plll'+ 

7 (2,0,0) 

-1 

1 

3 

4 

3 

0 

1 

3 






plll,+ 

7 (2,1,0) 

1 

0 

2 

3 

2 

3 

1 

3 

1 

3 





plll,+ 

7 (2,1,1) 

-1 

1 

3 

2 

3 

0 

1 

3 

2 

3 

4 

9 




7 (2,2,0) 

1 

2 

3 

10 

3 

4 

3 

2 

3 

2 

3 

8 

9 

4 

3 



T> lu t + 

7 (2,2,1) 

-1 

1 

3 

-2 

4 

3 

2 

3 

0 

0 

4 

3 

8 

9 


P 111 ^ 

(2,2,2) 

-1 

-1 

-6 

16 

3 

1 

2 

4 

-4 

8 

3 

16 

9 


Table 4. The coefficients of the polynomials P\ k \' k2 ks ^ with k\ < 2. 



1 

*i 

*2 

*3 

*i 

*1*2 

*1*3 

*2 


*3 

7 (0,0,0) 

1 










pill' - 
7 (1,0,0) 

-1 

1 









pill’ - 
7 (1,1,0) 

3 

-1 

2 








7 (1,1,1) 

-3 

2 

-2 

4 

3 







pill' - 

7 (2,0,0) 

-3 

-1 

-2 

0 

1 






pm- - 
7 (2,1,0) 

1 

1 

0 

4 

3 

-1 

2 





7 (2,1,1) 

-3 

-2 

-2 

0 

2 

-2 

4 

3 




pill’ - 
7 (2,2,0) 

8 

0 

12 

4 

3 

-2 

-2 

4 

3 

4 



pill' - 

7 (2,2,1) 

-6 

2 

-10 

4 

1 

2 

0 

-4 

8 

3 


P 111 ’- 

(2,2,2) 

-6 

-2 

-10 

16 

3 

3 

0 

16 

3 

-4 

8 

3 

16 

9 


Table 5. The coefficients of the polynomials P( k ^ k2 k3 ) with fci < 2. 


4.2. Continuous orthogonality. 

Continuous orthogonality of antisymmetric and symmetric cosine functions within each family is detailed 
in Section |T2| Our goal is to reformulate the orthogonality relations (18) - (21) after the change of variables 
(xi ,..., x n ) to polynomial variables (Xi, ..., X n ) and obtain the orthogonality relations for polynomials T k ± 
and In order to determine the corresponding weight functions in the integrals defining continuous 

orthogonality, we calculate the Jacobian of the change of variables (Xi(xi, .. 
variables (xi,..., x n ). 


X n (x 1 ,...,x n )) to 
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Proposition 4.2. The determinant of the Jacobian matrix J(x i,..., x n ) = det d( x i,...,x n ) for the coordinate 

v 7 0 \X l,---,x n ) 

change from (X 1 ,..., X n ) to (sq,..., x n ) is given by 


J(x\, . . . , X n ) = ( — 1) n< ” 2+1> C 


n(n-l) 


Tr n JJ(n — i)\i\ sin (lj2 ^ n) ( Xl ,... ,x n ), 


(73) 


\i= 1 


where 


n 

= 2 ” (Tl— si n (7ra:i) JJ sin (|(o: i + ^)) sin (|(a:j - xj)') . (74) 


i=l 




Proof. We prove the formulas (73), (74) by direct calculation. By the definition of the variables Xj via 
symmetric cosine functions, we have 

n 

x j = ( n ~j)'f- ^2 cos (nx^)... cos (nXij ). 

1 h<-<ij J 

Thus, the partial derivatives of Xj with respect to Xk are given by 

dX t 


j = 1 


3> 2: 


dx k 
d_Xj 

dx k 


= — 1)! 1! sin (Trx k ), 




-n(n — j)\j \sin (-KX k ) 


T. COs(wXi 1 )...COs(7TX ij _ 1 ) 
f — i = l 1 
\\ J 


— COS (7TXfe) 


y, cos (nx^).. .cos (nx i:i _ 2 ) 

il<’~<i 2 j 2 } 


■ COS ( 


; (nxk) — cos (nxk) cos ( 7ra ' i i) — cos ( 7TX k)J ■ ■ ■ 
Therefore, using the properties of determinants, we rewrite the Jacobian as 


J(x i, ...,x n )= det 


/ —7r(n—l)!l! sin(7rtci) 7r(n—2)!2! sin(7rtci) cos(7rtci) ... ( — l) n 7r0!n! sin(7rtci) cos n (vrtci) 
-7r(n—1)!1! sin(7TfC2) 7r(n—2)!2! sin(7rtC2) cos(7TfC2) ••• ( — l) n 7r0!n! sin(7r^2) cos n_1 ( 77 : 22 ) 


\ —7r(n—l)!l! sin(7rtc n ) 7r(n—2)!2! sin(7rtc n ) cos(7Tfc n ) ... ( — l) n 7r0!n! sin(7rtc n ) cos n 1 ('KX n ) , 


Comparing this form of the Jacobian with (73), we observe that 

/ sin(7rtci) sin( 77 : 21 ) cos( 77 : 21 ) ... sin(7ra:i) cos n_1 ( 77 : 21 ) 


sm 


(1,2,...,n) 


n(n-l) 

(xi, ... ,x n ) = 2 2 det 


sin(7TfC2) sin(7TfC2) COs(7TfC2) ••• sin(7TfC2) COS n ( 77 : 22 ) 
\sin(7rtc n ) sin(7rtc n ) cos(7Ttc n ) ... sin(7ra: n ) cos n_1 (7rtc n ) > 


Therefore, the equality (74) follows directly from the proof of Proposition 2.1 


□ 


One verifies that sin^ 2 ^ • sin^ 2 n ^ decomposes into a sum of symmetric cosine functions labeled by 
integer parameters 


Sim, d \ * sin,. <-) \ — 

(l,2,...,n) (1,2,. ..,n) o 


l^sgn(a) y (-1) 


a(ai,...,a n ) cog + 


creS n 


f CLi=± 1 \ 

\ i=l,...,n J 


(l+aicr(l),...,n+a r! ,cr(n)) ’ 


where a(ai,..., a n ) is the number of positive numbers among ai,..., a n . Therefore due to Proposition 4.1 
the square of the Jacobian J 2 is expressed as a polynomial p /,+ in X± t ♦.., X n of degree 2n, 

J(x) 2 =p I ’+(X 1 (x),...,X n (x)). 

The polynomial p I,+ , being the square of the Jacobian J, satisfies 

(*),..., *„(*)) >0. (75) 

Moreover, note that due to the equality ( |74| the Jacobian J does not vanish for any point in the interior of 
FiSff 1 ) and therefore 

p I,+ (Xi(x), ..., X n (x)) > 0 , xeF(Sf)°. ( 76 ) 
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The domain is the defined as the transformed domain F(*S^ ff ), i.e., 

d(Sf) = {(*!(*),.. .,X n (x)) I X e F(Sf)} . 

Fig. [i] shows the domain ) for n = 3. Note that the relation ( |75| implies that 

p I ' + (x 1 ,...,x n )> o, (x 1 ,...,x n )e^(Sf). 

Thus, we use the polynomial p /,+ to rewrite the absolute value of the Jacobian J as a function J of AJ,..., X 
on £(S£ ff ), 

J(X 1 , ...,X n ) = (X 1 ,...,X n )e $(Sf), 

i.e., it holds that 


(77) 


\J(x)\ = yV'RAl (*),..., *«(*))• 


Similarly, it follows from (67) that the products cos pi -cos pi , cos+ -cos+, and cos p2 • cos p2 are expressed as 
polynomials and J in ~ in Xi ,... ,X n such that 

J I ~(X 1 (x),...,X n (x )) = cos” (x)cos” (x), (78) 

J III ’ + (X 1 (x),... ,X n (x)) = cos+(x)cos+(x), (79) 

J III ~(X 1 (x),...,X n (x )) = cos~(x)cos~(x). (80) 

Example 4.2. For n = 3 are the explicit forms of the polynomial p /,+ determined by 

p I ’+(X 1 ,X 2 ,X 3 ) = 7r 6 (-8X 3 + X?X% - 12X§ + 12X 3 X 2 X 3 - A 3 X 3 ) ((3X 2 + 6) 2 - (X 3 + 3X r ) 2 ), 
and the polynomials J 111 ^, and J T11 '~ are given by 

J I ’-{X 1 ,X 2 ,X 3 ) = ] (-8X$+X 2 X 2 - 12Xl + 12X 3 X 2 X 3 - ^XfX 3 ) , 


J ni ’ + (X U X 2 ,X 3 ) = - (X 3 + 3X 2 + 3Xi + 6 ), 

J III ’~(X 1 ,X 2 ,X 3 ) = I(X 3 + 3X 2 + 3Xi + 6) (-8X 3 + X 2 X 2 - 12 X 2 + 12X!X 2 X 3 - ^X 3 X 3 J . 

In order to use the integration by substitution, we show that the transform 

p:xe F(Sf) -> (X 1 (x),.. .,X n (x)) e $(Sf) (81) 

is in fact one-to-one correspondence. Any continuous function / on F(S^) is expanded in terms of symmetric 
cosine functions |ll] and thus also in terms of polynomials V as 

/(x) = Y, B k V I k ' + {X 1 {x), ...,X n (x)). 


Assume that there are y,y G F(S^) such that y ^ y and (X\(y ),... ,X n (y)) = (X\ (5/),... ,X n (y)). Let us 
define a continuous function f{pc) = where the index i is chosen in such a way that yi ^ yi. Using the 
expansion in polynomials, we obtain that f(y) = f(y ), which is in contradiction with yi yi. Therefore, we 
conclude that the transform (p is injective. 

Due to the orthogonality relations ([l8|)-([2l]), Proposition |4.1| and integration by substitution (|8l|), we 

i_i I I j | j j j -j- >—• 

deduce the continuous orthogonality relations for polynomials V k ,V k ’ in the following statement. 

Proposition 4.3. Let X = (AJ, ..., X n ) and dX = dX 1 ... dX n . Each family of polynomials 

forms an orthogonal basis of the vector space of all polynomials in R[X] with the scalar product defined by 

the weighted integral of two polynomials /, g, 


L 


f{X)g(X)w(X)d,X , 


5(S£«) 

where w(X) denotes the weight function for each family of polynomials, 


w(X) = < 


•>!,+ 


w ,+ (X) = jpg for V k , 

^’■(U = A^T i for rl’~, 


///,+ 




W — J(X) 

(x) — 

— J(X) 




for Vi 11 *. 
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The continuous orthogonality relations are of the form 

J V I k ' + {X)V I k ; + (X)w I ’ + (X)dX =h k H k 5 kk ,, 

j v I k '-{x)v I k r{x)w I '-{x)dx =h k 5 kk ,, 

I r I k II ’ + (X)V I J I ’ + (X)w III ’+(X ) dX=2~ n H k 5 kk ,, 
J V I k II ’-(X)V I J I ’-(X)w in ^(X) dX=2~ n 5 kk ,. 


( 82 ) 


5. CUBATURE FORMULAS 


5.1. Sets of nodes. 

In Sections 3.2 and 3.4, the discrete orthogonality relations of the four symmetric cosine functions and 
the four antisymmetric cosine transforms over finite sets F^ ,=b ,..., F^ m,± are detailed. The orthogonality 

. Lv v,± s iven b y 


relations of the types I-IV, which are described in 


11 


are defined over the sets F^,..., ±■ N 


fU = <! ( A,..., A 


■ N 




{( 


ri 


N 

r\ 


r, 

’’TV 

1 

2 




pill,: 

■ N 


F IV >± - F 
r N ~ r N 


N 

M(A'A) .'") eD »} 

II, ± 


N 


») e £>t+l} 

(n,...,r n ) e D 


(83) 

Nj 5 (84) 

(85) 

( 86 ) 

The families of polynomials V and inherit these orthogonality relations which are crucial for the 

validity of cubature formulas. The sets of nodes arise from the finite sets ,..., F^ m,± as images of the 
transform p, given by (81). Denoting these images of the sets F^,..., F^ m,± as ... 

s {<p(s) I 8 <E F**} , t G {I,.. •, VIII} , 

we prove that the restriction of p on any set F^ is injective in the following Proposition. 

Proposition 5.1. 


^VIII,± ie 


(87) 


-it.i I 


te{I,...,VIII}. 


|UAT I — | r AT 

Proof. Suppose that there exists s £ F^ + such that s/s and (/?(s) = (/?(s), i.e., 

(X 1 ( 5 ),...,X n ( 5 )) = (X 1 (S),...,X n (5)). 

According to the expansion resulting from SMDCT I in |ll], any function / given on FP is written as a 
linear combination of a finite number of functions cos^ and therefore as a polynomial in Xi, ..., X n . Using 
this expansion in polynomials, we have for the function f(s) = (with i chosen in such a way that j-'sf) 

that f(s) = f(s) - which contradicts s^s. Thus, we obtain 
orthogonality relations, the proof for other grids is similar. 


5rh+ 

>N 


F I,+ 
r N 


. Using the corresponding discrete 

□ 


As an example, the points of the set for N = 5,10,15 in dimension 3 are depicted in Fig. [ij 

Due to the one-to-one correspondence given by p restricted to the grids F^, the symbols Hy, and £y 


are for p(s) = Y well defined by the relations 

PLy = H s , £y = £ s , 


£y = £ s - (88) 

has properties which yield to optimal Gaussian cubature formulas. 


5.2. Gaussian cubature formulas. 

Each family of polynomials 

Theorem 5.2. 

(1) For any N G N and any polynomial / of degree at most 27V — 1, the following cubature formulas are 
exact: 


F /0>'.+(y)<jy=(y)” £ 

A(ss ff ) 


WyVOO , 


(89) 
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Figure 8. The integration domain S f ff ) together with points of for TV = 5,10,15. 


[ f(Y)cu in ’+(Y)dY= (-T-) X! ^U/(W m ’ + T)- 
A(S£ ff ) \2N + 1 J f^ II + 

r N 


(90) 


(2) For any TV G N and any polynomial / of degree at most 2(AT — n) + 1, the following cubature formulas 
are exact: 


/. f(Y)cj I ’~(Y)dY= (-A £ f{Y)J^~{Y), 

h{s^) \ N J 

r TV 

/. f(Y)cj in ’~(Y)dY= (-7—) X! f(Y)J in ’-(Y)- 
JS(S-«) V 2iv + 1 / „,r^n.- 


(91) 

(92) 


Proof. We show that the formula (90) holds. Due to the linearity of the integrals and the sums, it is sufficient 

1. Such monomials are expressed as a product 
1. From Propositions |4.3| and |4.1[ 


to prove it only for monomials / of degree not exceeding 2TV - 
of two monomials - g\ of degree at most TV and g 2 of degree at most TV 

the polynomials g\ and g<i are rewritten as linear combinations of polynomials V k II,+ with ki < N and V l k , 
with k[ < TV — 1. Therefore, we only need to show the formula for those polynomials. Using the continuous 
orthogonality (82), we obtain 


L 


v, 


111,+, 


( Y)P I J I ’ + (Y)oj ni ’+(Y)dY = ^H k 5 kk , 


Jir t + 




Similarly, from equality (63) it follows that 
2 


27V + 1 


£ Hr 1 Vi II ’ + (X(s))Vif I ' + (X(a))(coa+(a))' 2 = fH k 5 kk , 


seF, 


VIII,+ 
N 


if ki < TV — 1 and k[ < TV — 1. Note that if k\ = TV, then T^ //,+ (F) = 0 for all points Y from ^ III,+ . 
This implies that we extend the discrete orthogonality relation for any k with k\ < TV and any k' with 
k[ < N — l. Consequently, the second cubature formula follows from the continuous and discrete orthogonality 
of polynomials We prove the other formulas similarly. Note that SMDCT II and AMDCT II needed 

to derive two of the results are found in 11 . □ 


Theorem 5.3. The cubature formulas (89)—(92) are optimal Gaussian cubature formulas. Moreover, it holds 
that 


• the orthogonal polynomials V k + with k\ = TV vanish for all points of the set ^ 

• the orthogonal polynomials V k ,+ with ki = N vanish for all points of the set $ N ,+ , 

• the orthogonal polynomials V k ~ with k\ = TV — n + 1 vanish for all points of the set 

• the orthogonal polynomials V k ~ with k% = N — n -\-1 vanish for all points of the set $ N 

Proof. The fact that the nodes are common zeros of the specific sets of orthogonal polynomials follows directly 
by substituting the grid points to th e d efini tion of the polynomials via symmetric and antisymmetric cosine 

we obtain that the number of polynomials V k + of degree N — l 


*n,+ 

>AT 


functions (66). By Propositions 


4.1 


and 
11,+ 


5.1 


is equal to the number of nodes in $ N 
other cubature formulas is similar. 


, and therefore the cubature formula (89) is Gaussian. The proof for 

□ 
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The Gaussian formulas (89)—(92) are special cases of formulas derived in a general setting as an (anti)symmetrization 
of orthogonal polynomials of one variable |2 . In particular, the variables ( 14 , ... ,u n ) from |2 are related to 
our variables (65) by the relations 

Xj 

U% (n — i)\i\ 


and our orthogonal polynomials (| 66 |) correspond to the measure of the form 

dy\ • • • dy n , -1 < y 1 < • • • < y n < 1 , 


IF 1 -ViT^ + yiY ( II (yi-Vj ) 2 

\l <i<j<n 


i= 1 


where 


uk=u k (yi,...,y n )= E Vh 

1 <i 1 <--<ik<n 


• Uik 


k = 1 ,..., n . 


The parameters a,/3 ,7 take for each class of the polynomials ( 66 ) the following values: 


V 


I,+ 


V 


7*- 


III,+ 
k 


7 III,~ 


a ~ P = 7 = , 

o 1 1 

a = P = --, 1 =- 1 

1 R 1 

a = 1 =--,P=-, 

1 R 1 

a=--,P= 1= -. 


The advantage of the current specialization is, besides more explicit formulation ready for application, the 
theoretical link between the polynomials and the (anti)symmetric cosine functions. This connection, similar to 
the connection of the classical Chebyshev polynomials of the first kind and the cosine function of one variable, 
may serve as a starting point for further investigation of the properties of the multivariate polynomials ( 66 ) 
such as their zeroes, discrete orthogonality, and Lebesgue constant. 


5.3. Other cubature formulas. 

Similarly as in Theorem |5.2[ one uses the remaining 12 discrete transforms to derive additional cubature 
formulas for the orthogonal polynomials. Note that these formulas are slightly less efficient than Gaussian 
cubature formulas. 


5.3.1. Formulas related to T^’ + . 

The transforms SMDCT I, V, and VI are used to derive additional cubature formulas for the orthogonal 
polynomials T^ ,+ . 


(1) 


For any N G N and any polynomial / of degree at most 2 N — 1, the following cubature formula is 
exact: 




E £YH?f(Y). 


Ye$ 


1,+ 

N 


(2) For any N E N, N > 2, and any polynomial / of degree at most 2(TV — 1), the following cubature 
formulas are exact: 


/_ f(Y)u I ’+(Y)dY= (—) E SrU^fiY), 
h(s^) \2M-1J + 

Y N 

/_ f(Y)u I ’+(Y)dY= (-T-) E Mrf(Y )• 


Yesff’t 


5.3.2. Formulas related to V T k ~. 

The transforms AMDCT I, V, and VI are used to derive additional cubature formulas for the orthogonal 
polynomials V k ~. 

(1) For any N G N, N > n, and any polynomial / of degree at most 2(TV — n) + 1, the following cubature 
formula is exact: 


L 


S(-S£ ff ) 


/( 5 > 




(y) dY 



E ^/(y)j 7 ’-(y). 
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(2) For any TV G N, TV > n, and any polynomial / of degree at most 2(TV — n), the following cubature 
formulas are exact: 


/_ f(Y)u I ^(Y)dY=(-^—) E 

Ms?) \ 2 M- 1 J 

Y N 

[ f (Y)u 3 i ’-(Y)dY=(-^—]) e £wmV’“en- 

Ms?) \ 2 N- 1 J 


5.3.3. Formulas related to . 


Vi 


The transforms SMDCT III, IV, and VII are used to derive additional cubature formulas for the polynomials 

///,+ 

k 

(1) For any TV £ N, TV > 2, and any polynomial / of degree at most 2(TV — 1), the following cubature 
formulas are exact: 


L 


3(S?) 


f{Y)t 




(Y)dY — ((.) E £ Y 'Hy 1 f{Y)j III ^{Y), 




/. /(y)w m -+(y)dy= (V E 

(2) For any TV G N, TV > 2, and any polynomial / of degree at most 2(AT — 1) — 1, the following cubature 
formula is exact: 

/ /(iV / 7 ’+(y)dy= (VT-A E £ Y n- l f{Y)j IH ^{Y). 

h(s^) \2N-IJ 


5.3.4. Formulas related to Vj/ 1 ' . 


V x 


The transforms AMDCT III, IV, and VII are used to derive additional cubature formulas for the polynomials 
iii- 


(1) For any TV G N, TV > n, and any polynomial / of degree at most 2(TV — n), the following cubature 
formulas are exact: 

/_ f(Y)u ni ’-(Y)dY= (V) E £Yf(Y)J IH ’-(Y), 




L 


f(Y)uj in '-(Y) dY = 


9(s?) 


N 


E f(Y)J in ’-(Y). 




(2) For any TV G N, TV > n, and any polynomial / of degree at most 2(TV — n) — 1, the following cubature 
formula is exact: 


L 


S(S* ff ) 


f(Y)u ni ’~(Y)dY 


2N-1 


E £ Y f{Y)J IH ’-{Y). 




6. Concluding remarks 

• In this paper, only the antisymmetric and symmetric generalizations of cosine functions are investi¬ 
gated. Similarly to the cosine functions, generalizations of the common sine functions of one variable 
are defined in 11 and these generalizations have several remarkable discretization properties. The 
discrete antisymmetric and symmetric sine transforms of types I-IV are developed in [8] for the two- 
dimensional case only and the discrete multivariate sine transforms of type I are found in [11 . The 
generalization of the antisymmetric and symmetric discrete sine transforms, analogous to the discrete 
multivariate sine transforms of types II-VIII 3], have not yet been described. 

• The 16 DCTs, which are described on the grids ,..., F^ 111, , are straightforwardly translated 

into 16 transforms of the corresponding polynomials on the grids, 

~i,± ~viii,± 
on j • • • > on 


Polynomial interpolation formulas, similar to formulas (3.30) and (64), can obviously be formulated. 
The interpolation properties of these polynomial formulas as well as the convergence of the corre¬ 
sponding polynomials series poses an open problem. 
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The question of whether one can introduce multivariate Chebyshev-like polynomials of the second and 
fourth kinds in connection with the generalizations of the sine functions forms another open problem. 
The decomposition of products of two-variable antisymmetric and symmetric sine functions from [8] 
indicates the possibility of construction of such multivariate polynomials. Moreover, the continuous 
and discrete orthogonality of the antisymmetric and symmetric sine functions further indicates that 
the corresponding cubature formulas, useful in numerical analysis, can again be developed. 

Since the Weyl groups corresponding to the simple Lie algebras B n and C n are isomorphic to (Z/2Z) n xi 
S n 110], it is possible to show that the antisymmetric and symmetric generalizations of trigonometric 


functions are related to the orbit functions - C-, S'-, S s -, and S^-functions - studied in [T 2 HT 3 H 22 ] 


For 

example, the symmetric cosine functions coincide, up to a multiplication by constant, with C-functions 
and the antisymmetric cosine functions become, up to a multiplication by constant, S z -functions in 
the case B n: and S s -functions in the case C n . The grids considered in this paper, which allow the 
eight types of the transforms for each case, are, however, different from the grids on which the discrete 
calculus of the orbit functions is described. Thus, the collection of the resulting cubature formulas is 
richer and includes formulas of the Gaussian type. 

The exploration of the vast number of theoretical aspects as well as applications of the Chebyshev 
polynomials [5,25 is beyond the scope of this work. However, since the basic properties of these 
polynomials - such as the discrete and continuous orthogonality together with the cubature formulas - 
are replicated for the multivariate symmetric and antisymmetric cosine functions, one may expect that 
many other properties such as existence of FFT-based algorithms and the Clenshaw-Curtis quadrature 
technique 28 will find their corresponding multidimensional (anti)symmetric generalizations as well. 
The work presented in this paper may represent a starting point for these open problems and further 
research. 
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